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Abstract 

We review modified F{R) gravity as realistic candidate to describe the observable uni- 
verse expansion history. We show that recent cosmic acceleration, radiation/matter- 
dominated epoch and inflation could be realized in the framework of _F(_R)-gravity 
in the unified way. For some viable classes of _F(_R)-gravity, the Newton law is re- 
spected and there is no so-called matter instability (the very heavy positive mass 
for additional scalar degree of freedom is generated). The reconstruction program in 
modified gravity is also reviewed and it is demonstrated that any time-evolution of 
' the universe expansion could be realized in F(_R)-gravity. These models remain to be 

realistic also in the presence of non-minimal gravitational coupling with usual matter. 
It is shown that same model which passes local tests and predicts the unification of 
inflation with cosmic acceleration also describes dark matter thanks to presence of 
, additional scalar degree of freedom and chameleon mechanism. 

00 

^ ■ 1 Introduction 
p 

. Modified gravity suggests very natural answers to resolution of several fundamental cosmological prob- 

' lems. For instance, the observable universe expansion history may be described by modified gravity. 

OO , Indeed, it gives very beautiful unification of the early-time inflation and late-time acceleration thanks to 

different role of gravitational terms relevant at small and at large curvature. Moreover, the coincidence 
problem may be solved in such theory simply by the universe expansion. Some models of modified gravity 
' are predicted by string/M-theory considerations. 

. From another side, dark matter may be described totally in terms of modified gravity. Moreover, 

' modified gravity may be useful in high energy physics, for instance, to solve the hierarchy or gravity- 

GUTs unification problems. Finally, modified gravity may pass the local tests and cosmological bounds. 
Usually the evolution of the universe can be described by the FRW equation: 

\h'-p. (1) 

Here the spatial part of the universe is assumed to be flat. We denote the Hubble rate by H, which is 
defined in terms of the scale factor a by 

H^'-. (2) 

a 

In III), p expresses the energy density of the usual matter, dark matter, and dark energy. The dark energy 
could be cosmological constant and/or a matter with 'equation of state (EoS)' parameter w, which is less 
than —1/3 and is defined by 

w^^. (3) 
P 

Instead of including unknown exotic matter or energy, one may consider the modification of gravity, 
which corresponds to the change of the l.h.s. in ([T|). 
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An example of such modified gravity pretending to describe dark energy could be the scalar-Einstein- 
Gauss-Bonnet gravity [T], whose action is given by 



(4) 



Here Q is Gauss-Bonnet invariant: 

g = R^- AR^.R^^" + R^^.p^R'^^P" . (5) 

Another example is so-called i^(i?)-gravity (for a review, see [2]). In F(i?)-gravity models'^], the 
scalar curvature R in the Einstein-Hilbert action 

•Seh = J d'^Xy/^R , (6) 

is replaced by a proper function of the scalar curvature: 

SpiR) = J d^x^F{R) . (7) 

Recently, an interesting realistic theory has been proposed in [3], where F{R) is given by 

FiR)^-^{R + fHsiR)) , /H5(i?) ^ - ■ (8) 

In this model, R is large even in the present universe, and fHs{R) could be expanded by the inverse 
power series of R: 

Then there appears an effective cosmological constant Aoff as Aoff = m^ci/c2, which generates the 
accelerating expansion in the present universe 

In the HS-model ([5]), there occurs a flat spacetime solution, where i? = 0, since the following condition 
is satisfied: 

lim fHs{R) = . (10) 

An interesting point in the HS model is that several cosmological conditions could be satisfied. 

In the next section, we review on the general properties of _F(i?)-gravity. After some versions of F{R)- 
gravity were proposed as a model of the dark energy, there were indicated several problems/viability 
criteria, which we review in Section O It is shown how the critique of modified gravity may be removed 
for realistic models. In Section [4l we propose models [4j and ^5], which unify the early-time inflation 
and the recent cosmic acceleration and pass several cosmological constraints. Reconstruction program 
for _F'(i?)-gravity is reviewed in Section [Sj The partial reconstruction scenario is proposed. Section six 
is devoted to the description of dark matter in terms of viable modified gravity where composite scalar 
particle from F{R) gravity plays the role of dark particle. Non-minimal modified gravity is discussed in 
section seven. Some summary and outlook is given in the last section. 



2 General properties of -gravity 

In this section, the general properties of the -F'(-R)-gravity are reviewed. For general _F(i?)- gravity, one can 
define an effective equation of state (EoS) parameter. The FRW equations in Einstein gravity coupled 
with perfect fluid are: 

P^^H\ p = -1^[3H' + 2h) . (11) 
For modified gravities, one may define an effective EoS parameter as follows: 

"'cff = -1 - 1^ . (12) 
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The equation of motion for modified gravity is given by 

^g^.FiR) - R^^F'iR) - g^,.UF'{R) + V^V,F'(i?) = -yTf^)^, . (13) 
By assuming spatially flat FRW universe, 

ds'' ^ -dt' + a{tf {dx'f , (14) 

i=l,2,3 

the FRW-like equation follows: 

= + 3 + F'{R) - 18 {ah'^H + HH^ F"{R) + p^m) (15) 

There may be several (often exact) solutions of Without any matter, assuming that the Ricci 

tensor could be covariantly constant, that is, R^u cx g^^, Ea. (|13p reduces to the algebraic equation: 

= F{R) ~ 2RF{R) . (16) 

If Eg. p^ has a solution, the Schwarzschild (or Kerr) - (anti-)de Sitter space is an exact vacuum solution 
(seefB and rcfs. therein). 

When F{R) behaves as F{R) (x i?™ and there is no matter, there appears the following solution: 

(m-l)(2m-l) 

H^- , (17) 

which gives the following effective EoS parameter: 

6m^ — 7m — 1 
3(m — l)(2m — 1) 

When F{R) (x i?™ again but if the matter with a constant EoS parameter w is included, one may get 
the following solution: 



^^H^^ (19) 



and the effective EoS parameter is given by 



w + 1 , , 

Weff = -1 + . 20 

m 

This shows that modified gravity may describe early/late-time universe acceleration. 

3 Problems with F(i?)-gravity 

Immediately, after the i^(i?)-models were proposed as models of the dark energy, there appeared several 
works [71 [5] (and more recently in [21 HH]) criticizing such theories. 

First of all, we comment on the claim in [71. Note that one can rewrite F(i?)-gravity in the scalar- 
tensor form. By introducing the auxiliary field A, we rewrite the action ([7|) of the -F(i?)-gravity in the 
following form: 

j d-^^^i^'i^) iR-A)+ F{A)} . (21) 

By the variation over A, one obtains A — R. Substituting A — R into the action (|2ip . one can reproduce 
the action in ([7|). Furthermore, we rescale the metric in the following way (conformal transformation): 

ff^.. ^e"g^, , a = -hiF\A). (22) 
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Hence, the Einstein frame action is obtained: 



Via) = e-5(e--)-e-/(5(e-^)) = ^-|^ (23) 

Here g {e~°') is given by solving the equation cr = — In (1 + f'{A)) = lnF'(j4) as ^4 = 5 (e""^). Due to the 
scale transformation (P^ . there appears a coupling of the scalar field a with usual matter. The mass of 
a is given by 



2 da^ 2 \F'{A) {F'{A)y F"{A) 

Unless Too- is very large, there appears the large correction to the Newton law. Naively, one expects the 
order of the mass rria could be that of the Hubble rate, that is, rria ^ H ^ 10^'^^ eV, which is very light 
and the correction could be very large, which is the claim in [7]. 

We should note, however, that the mass ma depends on the detailed form of F{R) in general [TT] . 
Moreover, the mass rUa depends on the curvature. The curvature on the earth i?earth is much larger than 
the average curvature i?soiar in the solar system and i?soiar is also much larger than the average curvature 
in the unverse, whose order is given by the square of the Hubble rate , that is, i?earth ^ -Rsoiar ^ H-^- 
Then if the mass becomes large when the curvature is large, the correction to the Newton law could be 
small. Such a mechanism is called the Chameleon mechanism and proposed for the scalar-tensor theory 
in [12 . In fact, the HS model [3j has this property and the correction to the Newton law can be very 
small on the earth or in the solar system. In the HS model, the mass is given by (see also jl3j ) 

~ 2n(n + l)ci (^j • ^^^^ 

Here the order of the mass-dimensional parameter could be ~ 10"^^ eV^. Then in solar system, 
where R - lO-^^ eV^ the mass IS given by ml - 10-58+3" eV^ and in the air on the earth, where 
R ~ 10-5° eV^ ml ~ 10-36+14" eV^ The order of the radius of the earth is lO^m ~ {IQ-^'^ eY)'\ 
Therefore the scalar field a could be heavy enough if n ^ 1 and the correction to the Newton law is not 
observed being extremely small. On the other hand, in the air on the earth, if we choose n = 10, for 
example, one gets the mass is extremely large: 

- 10^3 GeV - lO^^ x Mpianck • (26) 

Here Mpianck is the Planck mass. Hence, the Newton law correction should be extremely small. 

Let us discuss the matter instability proposed in [S], which may appear when the energy density or 
the curvature is large compared with the average one in the universe, as is the case inside of the planet. 
Multiplying g^^'^ with Eq. ([13]) , one obtains 



F('HR) ^ „^p„ , F'{R)R 2F{R) _ r. 



Here T is the trace of the matter energy-momentum tensor: T = T^,^)^- We also denote d"F{R)/dR" 

by Let us now consider the perturbation from the solution of the Einstein gravity. We denote 

the scalar curvature solution given by the matter density in the Einstein gravity by Rh ^ {k^ /2)p > 
and separate the scalar curvature R into the sum of Rb and the perturbed part Rp as R — Rt + Rp 
{\Rp\ <C l-Rfcl)- Then Eq.(p7)) leads to the perturbed equation: 

^^^^'^ ' □i?p + ^S^^P^bV'Rp + U{Rb)Rp . (28) 



3F(2)(i?b) 3F(^){Rb) Fi-^){Rb 



4 



Here U{Rh) is given by 



iFi^){Rbf 3F(2)(i?b) 3F(2)(ii'b)2 3^^(2)(i?6)2 ^ > 

It is convenient to consider the case that i?b and Rp are uniform, that is, they do not depend on the 
spatial coordinate. Hence, the d'Alembertian can be replaced with the second derivative with respect to 
the time coordinate: □i?p —dfRp and Eq. pS)) has the following structure: 

= -dtRp + U{Rt)Rp + const. . (30) 

Then if U{Rb) > 0, Rp becomes exponentially large with time t: Rp ~ Q^f^O^^ and the system becomes 
unstable. In the l/i?-model [14j . since the order of mass parameter is 



lO^^sec ~ (lO-^^eV) ^ , (31) 



one finds 



C/(i?,) = -i?, + ^^4-(lO"''sec) 



6^** \gcm 

R,^[lQSecy^(^^\ (32) 
^ ' \gcm~-^/ 

Hence, the model is unstable and it would decay in 10^26 gg^, ^fQj- planet size). On the other hand, in 
l/R + i?2-model [H], we find 

U{R^) ^ ^ > . (33) 

Then the system could be unstable again but the decay time is ~ 1, 000 sec, that is, macroscopic. In HS 
model [3], U{Rb) is negative 13;: 

U{R^) - r^— f < . 34 

Cin[n + I) 

Therefore, there is no matter instability [13]. 

Let us discuss the critical claim against modified gravity in [9l [10] . As shown in ([T6| , as an exact 
solution, there appears de Sitter-Schwarzschild spacetime in i^(_R)-gravity. The claim in [9jll0j is that the 
solution does not match onto the stellar interior solution. Since it is difficult to construct explicit solution 
describing the stellar configuration even in the Einstein gravity, we now proceed in the following way: 
First, we separate F{R) into the sum of the Einstein-Hilbert part and other part as F{R) = R + f{R). 
Then Eq.([T3]) has the following form: 

1 1 ^2 

= -\9t.. {f{R) + A) + i?^./'(i?) + gpu^nR) - VpV,f'{R) . (35) 

Here —A is the value of f{R) in the present universe, that is, A is the effective cosmological constant: 
A — — /(i?o)- We now treat the r.h.s. in ([55]) as a perturbation. Then the last two derivative terms in 
([55)1 could be dangerous since there could be jump in the value of the scalar curvature R on the surface of 
stellar configuration. Of course, the density on the surface could change in a finite width A as in Figure 
[1] and the derivatives should be finite and the magnitude could be given by 

d^^^. (36) 
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Figure 1: Typical behavior of R and p near the surface of the stellar configuration. 



One now assumes the order of the derivative could be the order of the Compton length of proton: 

- mp - 1 GeV - 10^ eV (37) 

Here rup is the mass of proton. It is also assumed 

i?- i?e lO^-^^eV^ , (38) 

that is, the order of the scalar curvature R is given by the order of it inside the earth. 
In case of the 1 /R model ^ , one gets 

2 4 

in 11 

Df{R) ^ V^V,f{R) ^ ^ 10-20 eV^ » Re . (39) 

Then the perturbative part could be much larger than unperturbative part in (|35p . say, R ^ Re 
lO^^'^cV^. Therefore, the perturbative expansion could be inconsistent. 
In case of the model [3], however, we find 

nfiR) ^ W,W.f'{R) ^ ^ 4 ~ 10-3-"" eV^ . (40) 



Then if n > 2, we find Df'{R), V^Vi//'(i?) ^ Re and therefore the perturbative expansion could 
be consistent. This indicates that such modified gravity model may pass the above test. Thus, it is 
demonstrated that some versions of modified gravity may easily pass above tests. 

4 Unifying inflation and cosmic acceleration 

In this section, we consider an extension of the HS model ,3J to unify the early-time inflation and late-time 
acceleration, following proposals [H [S] . 

In order to construct such models, we impose the following conditions: 



• Condition that inflation occurs: 



lim f{R) = -A, . (41) 



Here is an effective early-time cosmological constant. 
Instead of (|4ip one may impose the following condition 



lim f{R) = ai?" . (42) 

R — *oo 
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f{R) 




Figure 2: The typical behavior of f{R) which satisfies the conditions (I4ip . (|44|) . and (j45|) . 




Figure 3: The typical behavior of f{R) which satisfies the conditions (|42p . (|44| . and (|45)) . 

Here m and a are positive constants. Then as shown in p^ . the scale factor a(t) evolves as 

a(t) cx t'^" , /io = XT^^ , (43) 
3(w + 1) 

and WcS = —1 + 2/3/io. Here w is the matter EoS parameter, which could correspond to dust or 
radiation. We assume m ^ 1 so that H/H^ ^ 1. 

• The condition that there is flat spacetime solution is given as 

/(O) = (44) 

• The condition that late-time acceleration occurs should be 

f{Ro) = -2Ro , f'{Ro) - . (45) 

Here Rq is the current curvature of the universe and we assume Rq > Rq. Due to the condition 
(|45|). f{R) becomes almost constant in the present universe and plays the role of the effective small 
cosmological constant: A; ^ — /(i?o) ~ 2Rq. 

The typical behavior of f{R) which satisfies the conditions ((4T|) . ((44|) . and (|45)) is given in Figure [2] 
and the behavior of f{R) satisfying (|42)) . (|44)) . and (|45)) is given in Figure [3l 

Some examples may be of interest. An example which satisfies the conditions (HT1) . (UH), and is 
given by the following action 

n \2n+l , r)2n+l 

/(^) ^ r ° . • (46) 



/o + A{(i?-i?o)'"+V<'+^} 
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Here n is a positive integer. The conditions and (US)) require 



-2 — = 2^0 , — = . (47) 
Jo + Jl^o -/i 



One can now investigate how the exit from the inflation could be reahzed in the model (j46[) . It is easier 
to consider this problem in the scalar-tensor form (Einstein frame) in (j23[) . In the inflationary epoch, 
when the curvature R = A is large, f{R) has the following form: 



1 /o 



/(^)--- + 72l^- (48) 



Hence, one gets 



and 



(2n+l)/o 



(49) 



2n + l 



h ft V(2"+l)/o 

Note that the scalar field a is dimensionless now. Let us check the condition for the slow roll, \ V' /V\ <C 1 
Since 

Via) , / f!a 



( V 

r) ^'\[2n+l)h) 



V{a 

if we start with tr ~ 1, one finds 

V'{a) fRo 



(51) 



Via) VA.' ' 

which is very small and the slow roll condition is satisfied. Thus, the value of the scalar field a increases 
very slowly and the scalar curvature R becomes smaller. When a becomes large enough and R becomes 
small enough, the inflation could stop. Another possibility to achieve the exit from the inflation is to add 
small non-local term to gravitational action. 

We now consider another example, where f{R) satisfies the conditions (|42|) . ([44]) . and (|45|) [5]: 

•'^ ^ H-7i?" ^ ' 

Here a, /?, and 7 are positive constants and n is a positive integer. When the curvature is large (i? — > 00), 
f{R) behaves as 

/(i?) -i?" . (54) 

7 

To achieve the exit from the inflation, more terms could be added in the action. Since the derivative of 
f{R) is given by 

f'iR) = ^— 2 — > 55 

(1 + 7^?") 

we find the curvature Rq in the present universe, which satisfies the condition f'{Ro) — 0, is given by 

l/n 



and 



mo) ^ -2Ro -^11+ + j ■ ^ ^ 
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Let us check if we can choose parameters to reproduce realistic cosmological evolution. As a working 
hypothesis, we assume (i^jot 3> 1, then 

i?o ^ (/3/a7)'^'" , /(^o) = -2i?o ^ -/3/7 (58) 

We also assume f{Ri) ~ {a/j)R''j ^ Ri. Here i?/ is the curvature in the inflationary epoch. As a result, 
one obtains 

a - 2i?oi?o 2", /3 - ARlRo'^"'R'}-\ 7 - 2RoRo^''R'}-\ (59) 
Hence, we can confirm the assumption f3"f /a ^ f if n > f as 

Thus, we presented modified gravity models which unify early-time inflation and late-time acceleration. 
One should stress that the above models ([46]) and (I53p satisfy the cosmological constraints/local tests in 
the same way as in the HS model 0]. 

5 Reconstruction of F(i?)-gravity 

In this section, it is shown how we can construct F{R) model realizing any given cosmology (including 
inflation, matter-dominated epoch, etc) using technique of ref.[TS]. The general _F(i?)-gravity action with 
general matter is given as: 

S ^ f d''x^{F{R)+C^,ucr} ■ (61) 



The action ((6T|) can be rewritten by using proper functions P{(j)) and Q{(j)) of a scalar field (p: 

S = j d^X^ {Pi^)R + Q{<P) + ^matter} ■ (62) 

Since the scalar field (j) has no kinetic term, one may regard (p as an auxiliary scalar field. By the variation 
over (f), we obtain 

O = F'(0)i? + Q'(</)) , (63) 

which could be solved with respect to (/> as = 4>{R)- By substituting (p — 4>{R) into the action (|62|) . we 
obtain the action of F(i?)-gravity where 

F{R) = P{cP{R))R + Q{<p{R)) . (64) 

By the variation of the action ([62|) with respect to g^,y, the equation of motion follows: 

- -\9^.u {P^R + Qm - R^.uP{cp) + V^V,P(<^) - 5/..V2p(<^) -f ^T^, (65) 
In FRW universe (dH), Eg. ([55)1 has the following form: 

= -QH'P{<t>)-Q{^)-&H'^''^^^'^^ 



dt 

= (4ij + 6i7^)p(0)+Q(0)-H2^!M+4F^^+P (66) 
V / dt dt 

By combining the two equations in (|66p and deleting Qicp), we obtain 

^ 2^!«1 - 2H'-mm + 4i^P(0) + , + p . (67) 
dt'' dt 

Since one can redefine (j) properly as </> = (t){(p), we may choose to be a time coordinate: (j) = t. Then 
assuming p, p could be given by the corresponding sum of matter with a constant EoS parameters Wi 
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and writing the scale factor a{t) as a = ooe^'*-' (oq : constant), we obtain the second rank differential 
equation: 

' ^ - '^'(^^^ + 45"(0)P(^) + E (1 + -0 p.oao-^^^+-^e-3(i+^.).(« . (68) 

If one can solve Eq. (f68|) . with respect to P{<j)), one can also find the form of Q{(j)) by using ([66|) as 

Q(0) = -6 {9'{cf>) f P{<t>) - ^9'{4>)^ + E P.oao '^^+"")e-3(i+-)^(« . (69) 

Thus, it follows that any given cosmology can be realized by some specific _F(i?)-gravity. 
We now consider the cases that (|68p can be solved exactly. A first example is given by 

= 50 + f . (70) 

For simplicity, we neglect the contribution from matter. Then Eq. (j68|) gives 

d^P f , 9l\dP_2g, 

(j) ) d(j) 02 



The solution of (j7ip is given in terms of the Kummer functions or confluent hypergeometric functions: 

P = z'^FKia, r, , z'-^FK{a - 7 + 1, 2 - 7; z) (72) 

Here 



1 + ffl ± y/gl + lOffi + 1 
z = go(p , a = 



1 = l±- ^ , FK[a,r,z) ^ ^ — 7V—r- 73) 

2 Tn'+l --- 7 + n— 1 n\ 



Eg. ([70)1 gives the following Hubble rate: 
Then when t is small, H behaves as 



^ 7(7+ 1) •• •(7 + n- 1) n\ 

H = go + ^ . (74) 

H-j, (75) 
which corresponds to the universe with matter whose EoS parameter is given by 

w = -l + :^. (76) 

On the other hand, when t is large, we find 

H^go, (77) 

that is, the universe is asymptotically deSitter space. 

We now show how we could reconstruct a model unifying the early-time inflation with late-time 
acceleration. In principle, one may consider g((j)) satisfying the following conditions: 

• The condition for the inflation {t — cf) ^ 0): 

5"(0) = , (78) 

which shows that H{0) = ^'(0) is almost constant, which corresponds to the asymptotically deSitter 
space. 
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• The condition fot the late-time acceleration (at t = ^ to): 

g"(to)=0, (79) 
which corresponds to the asymptotically deSitter space again. 
An example could be 

9 (</>) = .90 + .gi ^ -^TTT- ■ (80) 

Here go, gi, and c are positive constants and n is positive integer greater than 1. Note that g'{4>) is a 
monotonically decreasing function of if < < to We also assume 

< go - ^ « .90 . (81) 
c 

One should note that g'{0) = go corresponds to the large Hubble rate in the inflationary epoch and 

g'(to) = go ^ " to the small Hubble rate in the present universe. It is very difficult to solve 
with (jSOp . so we expand g{(t>) for small (j). For simplicity, we consider the case that n — 2 and no matter 
presents. Then 

g{cb) ^ go -^cj^' + 0{cb* or gj) . (82) 



Hence, one gets 



P(0) ^ p, + p,eao^^^ 



3 .90 .9o 5oJ I .90 .9o 



50 



'O(gl) . (83) 



Here Pqi Pi, P2, and P3 are constants of integration. Using boundary conditions we can specify different 
modified gravities which unify the early-time inflation with late-time acceleration. 
We may consider another model: 

a + b(p'^ 

1 + C(f>'^ 

Here a, b, and c are positive constants satisfying the condition: 



9'(<P) - ^ • (84) 



- < a . (85) 
c 

In the early time (f> = t ^ 0, we find 

H{t) ^ g'{t) = a + (6 - ac) + O (u^) . (86) 

Then we may identify a as a cosmological constant which generates the inflation. On the other hand, in 
the late time (j> — t ^ 00, the Hubble rate H is given by 

i^(^)=5'(^) = - + ^-|^ + o(r^) , (87) 

C C(p 

which tells that the effective cosmological constant generating accelerating expansion of the universe 
could be given by b/c. 

When 41 = t is small, by comparing ([86|) with ([82| . we may identify 

go , ac-b^ . 
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Then by using (j83p . we find that the corresponding P could be given by 



P(0) = P^ + Ple'''^ -(ac-b) 



P 









|e"^4 







p' 



-t^e-^-Pl 
90 



(89) 



Here Pq , P{ , Pj, and P/ are constants of integration. On the other hand, when = t is large, we find 



P(0) = P(f + Pi^e^^/^ + 



a — b/c 



AP^ / d0'e^ 




(90) 



Here P^ , P/, and P^ are constants of integration, again. 

The important element of above reconstruction scheme is that it may be applied partially. For instance, 
one can start from the known model which passes local tests and describes the late-time acceleration. 
After that, the reconstruction method may be applied only at very small times (inflationary universe) 
to modify such a theory partially. As a result, we get the modified gravity with necessary early-time 
behavior and (or) vice-versa. 



6 Dark Matter from F(i?)-gravity 

It is extremely interesting that dark matter could be explained in the framework of viable P(P)-gravity 
which was discussed in previous sections. 

The previous considerations of P(P)-gravity suggest that it may play the role of gravitational alter- 
native for dark energy. However, one can study P(P)-gravity as a model for dark matter. There have 
been proposed several scenarios to explain dark matter in the framework of P(i?)-gravity. In most of such 
approaches [TB], the MOND-like scenario or power-law gravity have been considered. In such scenarios, 
the field equations have been solved and the large-scale correction to the Newton law has been found and 
used as a source of dark matter. 

There was, however, an observation [T7] that the distribution of the matter is different from that of 
dark matter in a galaxy cluster. From this it has been believed that the dark matter can not be explained 
by the modification of the Newton law but dark matter should represent some (particles) matter. 

It is known that P(P)-gravity contains a particle mode called 'scalaron', which explicitly appears when 
one rewrites P(P)-gravity in the the scalar-tensor form (|23p . In the Einstein gravity, when we quantize 
the fluctuations over the background metric, we obtain graviton, which is massless tensor particle. In 
case of P(P)-gravity, when one quantizes the fluctuations of the scalar field in the background metric, one 
gets the massive scalar particles in the addition to the graviton. Since the scalar particles in P(P)-gravity 
are massive, the pressure could be negligible and the strength of the interaction between such the scalar 
particles and usual matter should be that of the gravitational interaction order and therefore very small. 
Hence, such scalar particle could be a natural candidate for dark matter. 

In the model [3] or our models (UHl) and ([55]) . the mass of the effective scalar field depends on the 
curvature or energy density, in accord with so-called Chameleon mechanism. As our models (I46p and 
([55)1 describe the early-time inflation as well as late-time acceleration, the 'scalaron' particles, that is, 
the scalar particles in P(P)-gravity, could be generated during the inflationary era. An interesting point 
is that the mass could change after the inflation due to Chameleon mechanism. Especially in the model 
(|46|). the mass decreases when the scalar curvature increases as shown in (j49|) . Hence, in the inflationary 
era, when the curvature is large, one may consider the model where is large. After the inflationary 
epoch, the scalar particles, generated by the inflation, could lose their mass. Since the mass corresponds 
to the energy, the difference between the mass in the inflationary epoch and that after the inflation could 
be radiated as energy and could be converted into the matter and the radiation. This indicates that the 
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reheating could be naturally realized in such model. Let the mass of a in the inflationary epoch be m/ 
and that after inflation be mA- Then for N particles, the radiated energy En may be estimated as 



E ^ {mi ~ tua) N , (91) 

which could be converted into radiation, baryons and anti-baryons (and leptons). It is believed that the 
number of early-time baryons and anti-baryons is 10^*^ times of the number of baryons in the present 
universe. Since the density of the dark matter is almost five times of the density of the baryonic matter, 
we find 

mi > 10^°mA . (92) 

In the solar system, one gets A = i? - 10"^^ cV^. Then if n > 10 12 and A, - lO^''^^'*, the order 
of the mass m^ is given by 

to2^ 10239-295^10" eV2 , (93) 

which is large enough so that a could be Cold (non-relativistic) Dark Matter. On the other hand, in 
l/i?-model, the corresponding mass is given by 

-Vi^,-^-10-^'eV2. (94) 

Here fj, is the parameter with dimension of mass and /i ~ 10~33eV. The mass mi/^ is very small 
and cannot be a Cold Dark Matter. The corresponding composite particles can be a Hot (relativistic) 
Dark Matter but Hot Dark Matter has been excluded due to difficulty to generate the universe structure 
formation. 

In the inflationary era, the spacetime is approximated by the de Sitter space: 



1=1,2,3 

Then the scalar particle a could be Fourier-transformed as 



J d^kaik,t)e-'^'' . (96) 
Hence, the number of the particles with k created during inflation is proportional to e'^'^ . Here 
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Then if 

^ > 4 , (98) 

sufficient number of the particles could be created. 

In the original /(_R)-frame ([7]), the scalar field a appears as composite state. The equation of motion 
in /(i?)-gravity contains fourth derivatives, which means the existence of the extra particle mode or 
composite state. In fact, the trace part of the equation of motion (|13p has the following Klein-Gordon 
equation-like form: 

3V2/'(i?) = R + 2f{R) - Rf'{R) - K^T . (99) 

The above trace equation can be interpreted as an equation of motion for the non trivial 'scalaron' f'{R). 
This means that the curvature itself propagates. In fact the scalar field a in the scalar-tensor form of the 
theory can be given by 'scalaron', which is the combination of the scalar curvature in the original frame: 

a = -ln(l + /'(i?,)) . (100) 

Note that the 'scalaron' is different mode from graviton, which is massless and tensor. 
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Eq. (H5)) shows that the mass, which depends on the value of the scalar field cr, is given by 



1+1 



-^-^U^ (101) 



2 /o/ 2n+l \/ ff V"+-^.g:43 
/2 \2n + 2j V(2n+l)/o, 

If the curvature becomes small, a becomes large and decreases. Then the scalar particles lose their 
masses after the inflation. The difference of the mass in the inflationary epoch and that after the inflation 
could be radiated as energy and can be converted into the matter and the radiation. 
By substituting the expression of a (^5]) into poip . one obtains 



,2 



2(2n+l)(n + l)/o 



(102) 



Note that A corresponds to the scalar curvature. Let denote the value of A in the inflationary epoch by 
Ai and that after the inflation by Aa- Then the condition ([92|) shows 



'^'-] >10i°. (103) 



niA \AaJ 

For the model with n ~ 2, the condition ([5^ or (|103p could be satisfied if Aj/Aa > 10'^, which seems to 
indicate that the reheating could be easily realized in such a model. 

Now we check if the condition ([98|) could be satisfied. Note Hq ~ A^. Eq. (ll02p also indicates that in 
the inflationary era, where A = R ^ Ai, the magnitude of the mass is given by 

A 2n+l 

which is large enough and the condition is satisfied. Here Eq. ([17|) is used. Thus, sufficient number 
of (T-particles could be created. 

Let us consider the rotational curve of galaxy. As we will see the shift of the rotational curve does not 
occur due to correction to the Newton law between visible matter (baryon or interstcUer gas) but due to 
invisible (dark) matter, and the Newton law itself is not modified. 

Let the temperature of the dark matter be T = l/fc/3 where k is the Boltzmann constant. First, we 
assume the mass of the scalar particle a is constant. As the total mass of dark matter is much larger 
than that of baryonic matter and radiation, we neglect the contributions from the baryonic matter and 
radiation just for simplicity. We now work in Newtonian approximation and the system is spherically 
symmetric. Let the gravitational potential, which can be formed by the sum of the dark matter particles, 
be V{r). Then the gravitational force is given by .F(r) = —nidV{r)/dr. If we denote the number density 
of the dark matter particles by ri(r), in the Newtonian approximation, by putting = SttG, one gets 

J^{r) = / W^(s)ds (105) 

and therefore V{r) is given by 

/r 7 />s 
^ J u^n{u)du . (106) 

If one assumes the number density n{r) of dark matter particles could obey the Boltzmann distribution, 
we find 

n(r) = iVoe-'^'"-^''') . (107) 

Here A'o is a constant, which can be determined by the normalization. Using l|106p and (|107|) and deleting 
n(r), the differential equation follows: 

(r^V'ir))' = 47rG'm^iVor-2e-'^"-^('') . (108) 
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An exact solution of the above equation is given by 

As the general solution for the non- linear differential equation (|108p is not known, we assume V{r) could 
be given by ()109|) . Then the rotational speed v of the stars in the galaxy could be determined by the 
balance of the gravitational force and the centrifugal force: 

77^ 9777 

m,- = -Tir) = m^V'ir) - . (110) 

r fjma-r 

Here 'm^, is the mass of a star. Hence, 

- ^ . (Ill) 

that is, V becomes a constant, which could be consistent with the observation. 

For the dark matter particles from /(i?)-gravity, the mass rrio- depends on the scalar curvature or 
the value of the background a as in (|10ip . The scalar curvature is determined by the energy density p 
(if pressure could be neglected as in usual baryonic matter and cold dark matter) and if we neglect the 
contribution from the baryonic matter, the energy density p is given by 

pir)^m^n{r). (112) 

Therefore it follows 

TO<T = riia- {p{r)) = rua (man{r)) , (113) 
which could be solved with respect to ma'- 

rria — rria- {n{r)) . (114) 
Furthermore by combining (|107l) and (|114p . one may solve nia- with respect to V{r) and Nq as 

ma^maiNo,V{r)) . (115) 

Then (|105p could be modified as 

.F(r) = - ^'""^ (^o,^(r)) r ^^^2^^ (^^^ (^^g) 

which gives, instead of (jlOSp . 

{r^V'ir))' = 47rGm^ (iVo, V{r)) Aror2g-/3m,(JVo,v(r))y(r) _ ^^^^^ 

Ea. (|117p is rather complicated but at least numerically solvable. 

For the model (j46p . if the curvature is large enough even around the galaxy, the mass nia- is given 
by (|102p . The scalar curvature A = R is proportional to the energy density (since the pressure could be 
neglected), A (x p, and the energy density p is given by (|112p . Then 

J_ r 2{n + l){2n + l)/o \ ^ 
I f? 



n[r) ^ { \ (m.(.))-i^ . (118) 



Using (|107p . one also gets 

^•''-)° ,. ™oM«'A'.)-«^(s^iiif^^ (lis.) 

(2n + 3)/3m„(r) mo [ fi ' 
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Here mo has mass dimension. By substituting (|119p into (|117p . it follows 




77io J m^{r)^ \ mo J ma{r)^ 

(l - 1„ i!W I 4 I + "g" + 'V' X^'r^ („..(.))«fa ^ (120) 

V mo y m„(r)2J 2 [ /f J 

It is very difficult to find the exact solution of ()120p . although one may solve ()120l) numerically. Then we 
now consider the region where <^ mo but ln(mcr/mo) is slow varying function of r, compared with 
the power of r. In the region, we may treat In {m^/mo) as a large negative constant: 

In (m<,/mo) - -C . (121) 

Then the following solution is obtained: 

2(2n + 3) 
J. \ 2n + 5 

ma{r) = mo ( — 

2 _ 4(2n + l)(2n + 9)C , ^ izi+l f 2(n + 1)(27^ + l)/o , 

= (2n + 5)/3 ^"^^"^ I Tl ( ■ 

Note that ro can be real for any positive n. Eq. pi9p shows that 

2(2,1 + 3) 

2(2n+ 1) 1 ( r \ r 



V{r) = \ — ^- — In-. (123) 

^ ' 2n + 5 (3mo \ro J ro ^ ' 

Note that the potential (jl23p is obtained by assuming the Newton law by summing up the Newton po- 
tentials coming from the /(-R)-dark matter particles ('scalaron') distributed around the galaxy. Eq. (|122p 
indicates that the condition m^ ^ mp requires r ^ ro- Then by using the equation for the balance of 
the gravitational force and the centrifugal force, as in piOp , we find 

2n + 3) 

.ocy , (124) 

which is monotonically increasing function of r and the behavior is different from that in (jllip . If there is 
only usual baryonic matter without any dark matter, the velocity is the decreasing function of r, if there 
is also usual dark matter, as shown in piip . the velocity is almost constant, if dark matter originates from 
/(i?)-gravity, as we consider here, there is a region where the velocity could be an increasing function 
of r. Of course, one should be more careful as these are qualitative considerations. The condition 
m^r <C mo requires r ^ ro but in the region faraway from galaxy, the scalar curvature becomes small 
and the approximation p02p could be broken. Anyway if there appears a region where velocity is the 
increasing function of r, this might be a signal of /(i?)-gravity origin for dark matter. For more precise 
quantitative arguments, it is necessary to include the contribution from usual baryonic matter as well as 
to do numerical calculation. In any case, it seems very promising that composite particles from viable 
modified gravity which unifies inflation with late-time acceleration may play the role of dark matter. 

7 Non-minimal modified gravity 

In this section, we consider the theory with non-minimal gravitational coupling as an extension to the 
F(i?)-gravity. 

In [18 j . the non- minimal gravitational coupling of the scalar field </> was considered: 



'9 



(125) 
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Here A is a geometrical invariant like scalar curvature or Gauss-Bonnet invariant. In the FRW universe 
(|14p . by the variation over 0, one finds 



(126) 



Here g is a constant. By combining (|126p and the gravitational field equation, we delete the scalar field 
(p and obtain FRW-\ike equation: 



6 , 



Pi + P2 



Pi 
P2 



q 



h 



f2-6f.^iH + H^) 



7H^) + 6H 



h 



In case /2 = and /i = fns in 



6i//^'i? . (127) 
|), if the curvature R is small, the following solution is obtained: 

a cx , (128) 

which gives the effective EoS parameter WcS in (fT2|) as 

1 ~ n 



WcS 



1 



(129) 



Then the late-time acceleration of the universe Wcs < —1/3 occurs when n > 2. Furthermore, the 
effective phantom regime < — 1 appears when n < — 1. Within such theory for proper choice of 
gravitational part (like in previous sections) one achieves the unification of early-time inflation with 
late-time acceleration. Moreover, there is no extra problems to pass local tests. 

One may also consider the theory with non- minimal coupling with electro-magnetic field |19j : 



S 
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FjR) ijR) 

2k2 4 



(130) 



Here I{R) is a proper function of the scalar curvature R. Even if gravity is the Einstein one {F{R) 
if I{R) is chosen as 

liR) = 1 + fHs{R) , (131) 



ifHsiR) is given in 



I), the power law inflation occurs, where 
a{t) cx , 



(132) 



when R/vT? ^ 1. Note that the large-scale magnetic fields are generated due to the breaking of the 
conformal invariance of the electromagnetic field through its non-minimal gravitational coupling. The 
main conclusion is that adding the non-minimal coupling with gravity, the qualitative results of previous 
sections remain valid. On the same time, new bounds for non-minimal gravitational couplings only 
restrict such couplings to tend to constant (or, sometimes, to zero) at current epoch. The important 
point of non-minimal modified gravity is that as purely gravitational part any F(R) gravity[51 1201 121j 
may be discussed in cosmological context. The only condition is that it should be realistic, i.e. to pass 
the local tests and cosmological bounds. 



8 Discussion 

In summary, we reviewed -F(i?)-gravity and demonstrated that some versions of such theory are viable 
gravitational candidates for unification of early-time infiation and late-time cosmic acceleration. It is 
explicitly shown that the known critical arguments against such theory do not work for those models. In 
other words, the modified gravity under consideration may pass the local tests (Newton law is respected, 
the very heavy positive mass for additional scalar degree of freedom is generated). The reconstruction 
of modified F{R) gravity is considered. It is demonstrated that such theory may be reconstructed for 
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any given cosmology. Moreover, the partial reconstruction (at early universe) may be done for modified 
gravity which complies with local tests and dark energy bounds. This leads to some freedom in the choice 
of modified gravity for the unification of given inflationary era compatible with astrophysical bounds and 
dark energy epoch. Moreover, non-minimal gravitational coupling with usual matter may be successfully 
included into above scheme. As a final very promising result it is shown that modified gravity under 
consideration may qualitatively well describe dark matter, using the composite scalar particle from F{R) 
theory and Chameleon scenario. 

Thus, modified gravity remains viable cosmological theory which is realistic alternative to standard 
Einstein gravity. Moreover, it suggests the universal gravitational unification of inflation, cosmic accel- 
eration and dark matter without the need to introduce any exotic matter. Moreover, it remains enough 
freedom in the formulation of such theory which is very positive fact, having in mind, coming soon precise 
observational data. 
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